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CN . Abstract 

^ • In a recent work [1] we consider the topological expansion for the non-mixed observables (including the free 

energy) for the formal Cauchy matrix model. The only restriction in [Tj was the fact that all the branch points 

have to be simple. This excludes a very interesting situation not encountered in the literature before, namely 

the case in which two branch points merge in such a way that no cycle is pinched. In this work we focus on this 

, situation and derive new formulas that apply to non-simple and non-singular branch-points. 
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1 Introduction 

The fat-graph interpretation of matrix models has attracted a lot of attention since it was first noticed (see 
[H [3J m [S] for a short selection of early important works) . The Feynman graph expansion of matrix models 
produces graphs with certain characteristics depending on the model (like orientability or colorability of the 
graph), and they are weighted in such a way that one can organize the expansion in terms of many of its 
components and properties. The graphs are naturaly organized by the Euler characteristics of the Riemann 
surface on which it can be embedded. The free energy (the logarithm of the partition function) then is a 
generating function for the counting numbers of connected graphs that can be embedded in a compact surface 
with a given genus. In fact any observable of a matrix model has an interpretation in terms of combinatorics 
of graphs, for example the one-point resolvent generates counting numbers of connected graphs that can be 
embedded in a compact surface of a given genus with one border to which a given number of edges of the graph 
are connected. 

The topological expansion is an immediate consequence of the Feynman graph interpretation of the matrix 
integral, and so it is inevitably connected with the counting of graphs. That is one of the reasons why the 
topological expansion of matrix models has been so much investigated. 

In the past years a lot of progress has been achieved (see [6j [3 El [9l [1] and references therein), and the full 
topological expansion for a considerably large amount of models has been found. The whole scheme has been 
translated to a framework without matrix models as well |10| giving it an identity of its own and an interest 
beyond matrix models. This evolution has produced noticeable results in fields like String Theory [TTJ |T2] and 
algebraic geometry [13l [HI [T5j [16] . 

The whole construction is based on an algebraic curve (called spectral curve whenever a matrix model is 
related to it) under the generic assumption that all branch-points are simple. This constraint is sufficient to 
ensure that the algebraic curve is not singular, a case where the whole topological expansion stops making sense. 
However, it is not a necessary constraint: non-simple branch-points can appear without making the algebraic 
curve singular. 

In this work wc find the first example of such a situation in the case of the Cauchy matrix model introduced 
in |17| and explored further in |18| . In [1] it was studied as a formal model. The Cauchy matrix model with 
partition function 



consists of two positive definite hermitian matrices with external potentials Vi and V2 and an attractive interac- 
tion with each other determined by the determinant in the denominator. The eigenvalue density distributions 
in encoded in the cuts of a three sheeted covering of CU {00} (see figure [TJ. In particular, the eigenvalue density 
distribution for the first matrix pi (x) is encoded in the branch cuts joining the sheets Xi and Xq along pieces 
of M+ while the eigenvalue distribution for the second matrix p2{x) is encoded in the branch cuts joining the 
sheets X2 and Xq along pieces of M_ . 

In [T] the authors developed the topological expansion for the Cauchy matrix model with the condition that 
neither of the eigenvalue density distributions met the hard edge a; = 0. This constraint allowed the authors to 
remain within the framework of topological expansion that has become standard until now. In fact, the only 
situation in which this standard framework does not work for the Cauchy Matrix Model is when the density 
distributions of both matrices meet the hard edge at the origin in which case the branch point x = is no longer 
simple. The resulting algebraic curve in not singular due to the fact that the merging branch points connect 
different pairs of sheets, in other words, no cycle is pinched due to this branch point coalescence. 




(1.1) 
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Figure 1: Three sheeted covering of CU{cxd} related to the Cauchy matrix model. In blue the cuts corresponding 
to the first matrix eigenvalue density distribution, in red the cuts corresponding to the one for the second matrix. 

The main result of this work is the extension of the Eynard-Orantin residue formulas for the topological 
expansion to the situation we have described where the spectral curve contains one non-singular branch-point 
with branching number Ub = 2, i.e. the simpler non-simple branch-point. The non-singularity of the algebraic 
curve is consequence of the fact that no cycle degenerates when the simple branch-points merge at the hard 
edge to form the multiple branch-point, in other words the algebraic curve does not experience any critical 
topological change in the formation of the multiple branch-point. In terms of the recursive residue formulas for 
the topological expansion, the addition of a new four-legged vertex together with the standard three-legged one 
is the most remarkable result. Graphically, the recursion takes the form 



p 




that we prove in section 14.5.21 and express graphicaly in 14.5.41 This recursion equation gives us a glimpse of a 
larger structure for higher order non-singular multiple branch-points that we leave for a future work. 
The paper will be organized as follows: 

• In section [5] we introduce the model and explain very basic facts that are already known in the literature. 
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• In section |3] we introduce the two master loop equations for this model. One of them contains all the 
information about the algebraic curve of the model. We explain the structure of the spectral curve and 
introduce the necessary definitions. We also study how the curve is modified when we move in the moduli 
space. Finally we extract from the loop equations the recursion equations that is the starting point of our 
solution. 

• In section 2] we solve the recursion relations found in the previous section by using residue formulas and 
the Abel differentials on the spectral curve. Then we find how the variations on the curve affect the 
observables of the model. 

• In section [5] the topological expansion of the free energy of the model we found by inverting the loop 
insertion operator. 

• Appendix |B] gives some extra information about the meaning and how to derive the loop equations found 
earlier in the paper. 

• Appendix [Cl summarizes the definitions and properties of all the algebro-geometric objects wc need in the 
topological expansion construction. 

• Appendix ID] contains the detailed account of some of the more tedious calculations. 



dP(Mi, A/2) = —dMidAh .jv ■ (2-2) 



2 The Cauchy Matrix Model 

The model of interest here is the Cauchy Matrix Model introduced in [T71 [IS] and whose topological expansion 
was found in [T] for all situations except a small subset. The partition function for this model is 

dMidMa ^ = e"^-^, (2.1) 

M, =Mj >o det (Ml +M2) 

where the domain of integration (as indicated) is the ensemble of N x N Hermitian positive definite matrices. 
Consider the matrices Alj as random variables with the probability measure 

1 p-f tr(yi(A/i)+y2(A/2)) 

2c det (Ml + M2 

The potentials Vi parametrize the measure associated to each matrix, and must satisfy a suitable growth condi- 
tion at infinity for the integral to be convergent; this condition is usually expressed as Vi{x)/ln{x) — > +00. 
However, we are interested in the formal expansion of the integral (the topological expansion) so this condition 
can be disregarded completely. 

We restrict this paper to the class of polynomial potentials of degree di 

di 

^^(^) > 4?>0, «-l,2, (2.3) 
fe=l 

although later we allow formal variations with respect to the infinite number of parameters t^j!^ . The logarithm 
of the partition function Z is called the free energy J- of the matrix model and the parameter T is called the 
total charge of the model. 

With the Harnad-Orlov formula [TO] and following (TB] we can write (|2.ip in eigenvalue representation 

Zc= l[dx.,dy., A^x,)A^y,) — , (2.4) 
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l<ij<JV 



where Xi (resp. yi) are the (positive) eigenvalues of Mi (resp. M2), and A(ai) = det ^ 

ni<j('^i — Gj) is the Vandermonde determinants. 

In the large limit the eigenvalues distribute along segments of the real positive line with eigenvalue 
distributions 

1 ^ 1 ^ 

pi(x) = Jim^-^5(x-x,) , p2{y)= \im^-Y,Hy-y^)■ (2-5) 

i=l i=l 

In the generic case the eigenvalue distributions may reach the hard edge at x = 0. This implies a generic behavior 
Pk{x) ^ x^i . When both eigenvalue distributions reach the hard edge, the attraction to the eigenvalues of the 
other matrix modifies this law to pk (x) ~ a;~ ^ . 

Note that we can easily rewrite equation (|2.4|) back into matrix representation with a different (but equivalent 
as long as we look at non-mixed correlators) interaction, namely 

r g-f tr(yi(Mi) + V2(Af2)) 

where now the determinant is acting on a N'^ x N'^ size matrix, and I represents the N x N identity matrix. 
Remark that the exponent TV in the determinant in (|2.1[) is not present in (|2.6p . 

We will consider here the formal perturbative expansion of the matrix integral around one local minimum 
generically constituted of a finite number of disconnected supports for the eigenvalue distributions of both Mi 
and M^. We don't extend here defining precisely the formal perturbative expansion of matrix models since it 
has already been defined in great generality elsewhere -see [?, ?] for more details-. For the sake of completeness 
we reproduce the main idea of the formal matrix model definitions in appendix [A] 

The goal of the this paper is to extend the results of [I] to the case where the supports of the eigenvalue 
densities for the two matrices meet the hard edge at the origin. 

3 Loop Equations, spectral curve and recurrence equations 

The two loop equations essential for the Cauchy matrix model have been derived in [1]. We have reproduced 
the main steps in appendix |BJ Following the indications in the appendix and in [T] the reader can reproduce 
these loop equations. 

3.1 Definitions 

In order to write the loop equations in a manageable form we need some definitions. Let us define the matrix 
average notation 

1 r p-f tr(Vi(Afi)+y2(M2)) 

= — / dMM (...) , (3.1) 



2cJau>o ' det(Mi (^I + I®M2)' 



^Note that we don't ask it to be the absolute minimum 
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We also define the one point resolvent functions, the t/-potentials and the F-functions, the main objects of the 
next section 



Wi{x 
Yiix 

W2{X 



T 
N 



1 



X - Ml ^ 
U[ix)~Wiix), Uiix) 
1 



2V(ix) - V^{~x) 



T 



(3.2) 



-V{{x) + 2Vi{-x) 



^X + M2, 

U^{x)~W2{x), U^{x) 
-Wi{x)-W2{x), 
~Yi{x) - Y2ix) ^ U'^{x) - Wq{x), 
where Uq is defined implicitely through the definitions of Yi , I2 , U'l and U2 ■ The one point resolvent functions 



(3.3) 



Y2{X 
Wo(.T 

Y^{x 



are related to the eigenvalue density distribution bj 

2TTipi{x) = lim lini (Wi{x - ie)) - (Wi{x + ie)) , 
2TTip2(—x) = - lim lim (W2(x — ie)) - (W2(x + ie) 



3.2 The loop equations 

The first of the two master loop equations (whose derivation is more explicit in appendix |B]) is quadratic in 
both Yi and Y2 and has the form 



(Yi{x))') + {{Y2{x))^) + {Yi{x)Y2{x)) = (i?(x-)) , 



(3.4) 



where {^R{x)j is an unknown meromorphic function of x with a pole of order 1 at most at a; = and a 
polynomial part, as can be seen from the explicit expression (see appendix IB|) 



R{x) = \{Vi{x)^ + Vii-x)^ ~ Vi{x)Vi{-x)) - Pi(x-) - P2{x) ^\{h+ h) 
= {U[{x))^ + (C/^(x))2 + U[{x)U'^{x) - Pi(x) - P2{x) -\{h+ h) , 



(3.5) 



where the polynomials (^Pi{x)J and (^P2{x)j as well as the constants \ and (^hj are defined in terms of 
matrix averages in appendix [Bl For future purposes we rewrite equation p.4p using the definition of Yo{x) 

2 

Y,(y,{x)')=2(r{x)). (3.6) 

A:=0 

The second of the two master loop equations is a cubic equation that both Yi , Y2 and Yq satisfy identically 

{y,{x)') {Rix)n{x)) ^ + i^) {W,{x)) = (Dix)) , for k = 0, 1,2, (3.7) 

where ^Z?(x)^ is a meromorphic function of x with a pole of order at most 2 at a; = and a polynomial part 
and can be expressed as 



D{x)) = ~U;,ix)U[{x)U!,{x) - U[(x) (P2ix)) - C/^(a;) (Pi{x)) + {Qi{x) + Q2{x)) 



(IKi) - U[{x) i^h) + U'2{x) i^h)) + ^ ( 



(3.8) 



''Note that for W2{x) is defined in such a way that the resolvent function's poles arc all on the negative axis and their residue 
is negative. This explains the signs in the definitions of p2- 
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where the polynomial (^Qi{x)j and \ Q2{x)j and the constants and \K2j are defined in terms of matrix 

averages in appendix[Bj Equation p.7p is essential to obtain the spectral curve that determines all the properties 
of the model. If we combine the cubic equation for all three we obtain 

2 

Y^lYr.{xf)=il^D{x)). (3.9) 

fc=0 

This equation together with equation p.6p is all we need to produce the recursion equations that allow us to 
compute the topological expansion. 

3.3 Spectral curve 

From equation p.7p we can extract the spectral curve that defines the leading order of (Yk{x)'^. Consider the 
leading terms in the topological expansion of the objects in (|3.7p 



n2 ■ 



(i?(x))=i?W(x) + 0(^^), (3.10) 
D{x))=D('\x) + o(^^y 



and extract the leading part of the equation p.7p 

yk{xf - R^°\x)yk{x) ^ d(''\x) fori = 0,1, 2. (3.11) 

The three functions yk{x) are the three solutions to the algebraic curve E{x,y) = 2/^ - R'-"\x)y - = 0. 
This algebraic curve is called the spectral curve of the matrix model. The algebraic curve can be parametrized 
by a Riemann Surface E. For every solution of the equation E{x,y) ~ there is a point p £ S, thus there is 
a parameterization x{p) and y{p) for the ensemble of solutions. Since there are three possible solutions of the 
algebraic curve for a given x, the curve can be represented as a three sheeted covering of C U {oo}. In other 
words, for a generic x, there are three and only three different points p,q,r € E with the same x-projection 
x{p) = x{q) ~ x{7') = X, and whose j/-projcctions are yk{x), k = 0,1,2. Reversing this argument we can 
construct the mcromorphic function on the algebraic curve 

V{p) = yk{x) , whenever p e A'^, (3.12) 

where is the x-based sheet of the algebraic curve defined by yk{x). From the definitions p.2|) we can 
see that y{j)) have cuts joining Xi and Xq along the positive real line (those coming from the eigenvalue 
distribution of Mi), and cuts joining X2 and Xq along the negative real line M_ like represented in figure 
[2 We calE Ap\ i ~ 0, . . . ,di — 1 the cuts corresponding to the eigenvalue distributions of Afi, and we call 
i = 0, . . . , d2 — 1, those corresponding to the eigenvalue distributions of M2 with dk < dk- We call a,|'^'' 
the branch-points and af'^ = x(q;-'^'') its x-projection, so that A-'''' = [a2i\a2i+i] support for the cuts. 

Denote by A['^\ i = 0, . . . ,dk — I, k = 1,2 small cycles on Xk containing A[''\ and cycles that go 
through Aq*^-* and Af'^ and whose x-projection encircles the branch-points aj*'', . . . ,a'^K The cycles A[''\b^'^'' 
-for i ~ 1, . . . ,dk — l, and k ~ 1,2- form a canonical homology basis on the algebraic curve with the intersection 
property Af M i5^'^ = S,,j6k.i (see fig.©. We are interested in a situation where cuts of both species reach the 



^Note that though the maximum genus is reached with di cuts A^^^ and d2 cuts a'^\ we allow a smaller number of cuts. 
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Figure 2: Three sheeted covering of C U {00} with a branch point at a; = with branching number 2. Red 
cycles are the A cycles, blue cycles are the B cycles. 

hard edge at 2; = like in fig. [21 It is clear from our definitions that the convergence of two simple branch-points 
and the corresponding formation of a branch-point with branching number 2 (connecting all three sheets) do 
not imply the degeneration of a cycle or a change in the genus of the Riemann surface g = di + d2 — 2. 

To specify one minimum for the topological expansion it is sufficient to fix the number of eigenvalues 
that belong to each disconnected component of the eigenvalue support. This is called the fixed filling fraction 
condition. The expansion itself cannot modify the minimum originally chosen, or equivalently the filling fractions 
do not receive corrections in N. The leading part of the fixed filling fraction condition can be written as 

/ y(p)da;(p) = 27rzef' , fc = 1, 2 , ^ = 1, . . . , 4 - 1, (3.I3) 

Ja\' 

and the rest of the expansion of the fixed filling fraction condition is 

/ w[''\p)=0 , fc = l,2,z = l,...,4-l, (3.14) 

where w['^'^ (p) -that we define in more detail later- is related to the h order term in the topological expansion of 
Wkix)y Ahhough i?(o)(a;) and D^°^x) are unknown a priori, choosing values for the filling fractions should 
amount to determine those two functional. 

For convenience we denote by p''' (i = 0,1,2) the three different points of the curve with the same x- 
projection as p, i.e. x(p^'-') — x{p) (obviously p = p'*^ for some value of i. We also denote the three points that 
project to a; = c« by ooi G Xi for j = 0, 1, 2. 

Consider also automorphisms defined locally around a branch-point a by 




(3.15) 



^We live it as an unproven statement here that the amount of unknown independent parameters in _R(") and DC'' are exaetly 
equal to di + 6,2 — 2 and that the eonditions 1 13.131 1 plus the faet that the genus \s g = d\ + d2 ^ 2 are enough information to 
determine all these unknowns. 



9 



that results from moving the point p around the branch point a in such a way that its x-projection turns 
counterclockwise exactly i times around a ~ x{a) and such that a;(i?^(p)) = x{p). It is easy to see that 
'&aip) = p if the branching number is 1, — p if the branching number is 2, and so on. We will drop the 

subindex a from now on unless the expression becomes ambiguous. 

3.3.1 Moduli of the model. 

The algebraic curve is parametrized by the data from the model 

(k) 

• The potential parameters tj ; 

• The total charge T; 

• The filling fractions e,|'^'' defined below in cq. (|3.22p . 

(k) (k)^ 

We denote a generic setting of these parameters as = {tj ,T,e\ }. We will be interested in computing 
variations (derivatives) with respect to these parameters. For the potential parameters, we introduce a gen- 
erating operator of such variations called the loop insertion operator in the matrix model literature. These 
operators are very important as they define the one point resolvent when applied to the free energy, as well as 
the connected n + 1-point resolvent when applied to the connected n-point resolvent. In our model we have two 
of them corresponding to the two sets of times , i = 1, 2 



y— 



dVi{x) ^ 9^(1) ' 

^^^^ 

dV2{x) ^ {-xy+'' Qt(^)' 

k k 

With these two operators we will define the global loop insertion operator acting on functions on the algebraic 
curve 

^ whenever p Cz Xi 

whenever p Cz Xq 
whenever p (z X2 



-whx)-dmT) whenever p e A-o (3.17) 



dV{p) 



dV2{x) 

We define as well the functiorlfl U'{p) on E 

U'{p) = U'kix{p)) whenever p e Xk. (3.19) 

The loop insertion operators p.l6p can be used to generate all the connected resolvent observables from the 
free energy T 



k=l 
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Remark that from the above definitions it follows by a direct computation that 



^ ;U'ip') = l , , , \ (3.18) 



dV{p) ' ' 3 (x{p) - x{p'))2 ■ 
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3.3.2 Fundamental differential of the algebraic curve 

The algebraic curve can be parametrized equally well by A4 or by the meromorphic functions y{p) and x{p). In 
fact we will see that the differential y{p)dx(j>) on E contains the information in M. Indeed,we can extract t[r''' 
for fc > from the behavior of y{p)dx{p) arounj^ 001.2 

yip)dx{p) ^ 2tW-(-l)4f) ^ 
p-^ooi (a;(p))*-' 3 ' 

y{p)dx{p) ^ -tj'^ + (-1)^24^^ 
3 



Res 



Also the filling fractions can be explicitly written as the A cycles of the fundamental differential y{p)dx{p) as 
we have already seen above 

^J^,^,yiP)Mp)=ef = i-ir^, (3.22) 



where ivj*^' is the number of eigenvalues belonging to the cut A''^\ The sign in p.22p appears due to the 
definition of W2 {x) . The total charge T can be represented in two different ways 



Res yip)dx{p) -T, 

(3.23) 

Res y{p)dx{p) — ~ T. 

p— >002 



The same expression evaluated at ooq reflects the condition X]i=o y(P 

Res y{p)dx(p) =0. (3.24) 

3.4 Variations of the spectral curve with respect to the moduU 

Consider the formulas from the previous section, it is clear that when we modify the moduli of the spectral 
curve, these changes must be reflected on the fundamental differential y(j))dx{p) . We are going to study how 
y(jj)dx{p) changes under the variations of the moduli. Since we have the constraint E{x,y) ~ 0, x{p),y{p) and 
p are not independent, we have to keep fixed one of the three while computing the variations. For convenience 
all variations are taken at x{p) fixed. We will generically call this variations dil. The variations will be written 
in terms of standard objects in algebraic geometry that we define in appendix [Cl 

3.4.1 Variation with respect to the filling fractions 

Not all filling fractions are independent due to the restrictions pointed out in the previous section. As with ^j*^' 

we consider as independent the for j = 1, . . . , dfe — 1. Since these are independent we have — fc-j- = Si jSk k' 

ae^ ' ' ' 

for i, j ^ 0. 

Let us compute the variation of y(p)dx(p). We have that 

y{p)dx{p) = -Tj-AU\p) - W{p))dx{p) = O {x{pr^) dx{p) (3.25) 



9f d,f 



^The behavior around the third infinity point p = ooq is not independent due to the relation y\{x) + y2{x) + ^3(2:) = 0. 



11 



when x{p) oo. Indeed the potential U'{p) docs not depend on the e's by hypothesis, and W{p) gives us the 
behavior at oo. On the other side we can compute the A cycles of the above differential 

'('fc) ' (3.26) 



'j,j'Ok,k' 



All together we have that —jp^y{p)dx{p) is a differential on the algebraic curve with no pole and with normalizing 

condition p.26p . this determines dfl (k) completely to be a basis of normalized first type Abelian differentials. 

j 

They can be written in terms of the Bergman Kernel B{p, q) 

dn^w = -l-y(p)dx{p) ^ 2Tnduf^ = / Bi.,p). (3.27) 

3.4.2 Variation with respect to the total charge T 

Imagine now that we change the total charge T, the variation of y{p)dx{p) will be 

r O {x{p)-'') dx{p) peXo 
^y{p)dx{p) = J il^) + O Hp)'^)) Mp) P^Xi (3.28) 
^ + 0{x{p)-^))dx{p) peX2 

Indeed by hypothesis the potentials do not depend on T, and the behavior at infinity of Wk{x) gives us the 
above behavior. This differential has poles only at cxji_2 of order 1 with residues +1, —1 respectively. Also the 
A cycles of this differential arc zero because the filling fractions do not depend on T cither. The conclusion is 
that _ 

d^T^p) - —y{p)dx{p) = d5ooi,oo.(p) = / B{.,p), (3.29) 

which is a normalized third type Abelian differential. 

3.4.3 Variation with respect to the potentials 

Consider now the potentials. The variations with respect to the parameters of the potentials look like 

dn^^rM EE -l^y{p)dx{p) = -^AU'{p) + O {x(p)-^))dx{p) 

ljx{py-^dx{p) + O {x{p)-^) for pCzXi 
-\jx{py-^dx{p) + O (.t(p)-2) for p(^X^ ' 



dr!^(.)(p) EE -^y{p)dx{p) = -^(C/'(p) + O {x{p)-'))dx{p) 



(3.30) 



, , \3^{py-^dx{p) + O (.t(p)-2) for p^X^ 

I -\3x{py~^^x{p)^0{x{p)-^) ioTp(^X2 ■ 

This behaviors at oofc plus the normalizing conditions /^(fc) dri^'^ = determine completely dfij*^'' as a combi- 
nation of second kind differentials 

d^lf\p) = \Rcsx{-yB{;p)~\ Res x{-yB{;p), 

■> 6 COi 6 002,000 

dnf\p) = {-ly-^ [-Resx{-yB{-,p) - - Res x{-yB{-,p) 
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Explicit vs implicit derivatives of y{p)dx{p). We claim that we can represent y(p)dx{p) as 



y{p)dxip) =2mY, E ^f^^t\p) 

k=l 3 = 1 



fc=i j 



(3.32) 



(k) (k)^ 

where ta is a collective notation for the whole set of parameters of our model A4 = {tj ,T,ej }. To see this, 
we provisionally denote by r](j)) the right-hand-side of (|3.32p : it is seen, from the definition of the various terms 
that ri{p) — y(p)dx{p) has no poles at cx)o, 001,002 and hence it is a holomorphic differential. On the other 
hand, from the definition of the filling fractions ej'^-' it also appears that it has vanishing y^-cycles. Standard 
theorems then guarantee then that ri{p) — y{p)dx{p) = 0. Therefore (|3.32p yields an alternative representation 
of the fundamental differential y(j))dx(j)). 

On the face of it appears that dt^y{p)dx{p) is the same as the coefficient of the corresponding "time" ta 
in the r.h.s. of (|3.32p . namely, the "explicit" derivative of the r.h.s. of (|3.32p . Note that this is not a trivial 
statement, since the spectral curve and all the differentials appearing as coefficients of the times in (|3.32p do 
depend (implicitly) on the times as well. 

Note that all variations dVlt^ are expressible as an integral operator on the Bergman kernel B{p,q) 



drtt^ = -^y{p)dx{p) = / Bip, ■)AtJ-) 

Ota Jct„ 



(3.33) 



where Ct^ is the integration path associated with the variation with respect to ta and Af^ (p) is a mcromorphic 
function on the curve. For example, looking at the formula above, we have Cji = fCooi — ^ (Cooo +^002)7 and 

Aj(i)(p) = x{py , where Coo^ is a small cycle around oofc. We thus introduce the integral operators associated 

j 

to the various derivatives 



dta 




Ja(n) 






pRes ^ 


Res ) x^D 




y 3 001 3 


000,002 J 


<>° 




(-^Res + 
\ 3 002 


- Res ] i^xya 

3 000,001 J 


34"° 




/ □ 




An 

dT 




} ° 

J0D2 





(3.34) 



where the regularized integrals J( arc defined (if necessary) by subtraction of the singular part in the variable 
x{p) (see [20]) With this notation we have -compactly- 

yip)dxip) = iaJa{B{-,p)). (3.35) 



3.5 Topological Expansion of the Loop Equations 

The equations we need to work with are equations p.6p , p.9p and the equation 

2 

E(n.(x)) = o. 

fc=0 



(3.36) 
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In terms of connected components we have 

2 



k=0 

2 / 2 y2 



E (^fc(^)) + {Wkixf) = 2 {Rix)) , (3.37) 

k=0 ^ 

j2((n{x)y + 3^{n{x)){w,{xr) ^^^(W,{xf) )=2,{d{x)] 

k=0 ^ c cj 

To proceed we need to define some more objects: write the topological expansion of our main observable^ 

oo / „ \ 2/i 



r,(x))=y,(x)-^f|) W\''\x\ 

h=l ^ ^ 
°o / rp\2h 



/i=0 



Ww.M) =Eh^ <';^,,„(xi,...,x„), (3.38) 

fe=l / c '1=0 ^ ^ 



/i=0 

rp\2h 



^(-) =E 4^ ^''(-)- 



/i=0 ^ 

From the /i components of these topological expansions we construct the following multi-differentials on the 
algebraic curve: 

y{p)dx{p) = y.i{x{p))dx{p) Vp e X^, 

j[''\p) = W^^'''{x{p))dx{p) ypeXi, h>0, 



ah).,_„Ah), ...... ^ _ „ (3-39) 



t«2^^(Pi,P2) = Wi'^l^{x{pi),x{p2))dx(jpi)dx{p2) Vpi G Xi^ andp2 G X^^, 
are the one and two-point resolvent differentials, while in general we define 

n 

W^J:\pi,...,Pn)=Wl';;l^^^{x{pi),...,x{pn))\{dx{pk) Vpfc G A',, , (fc=l,...,n) (3.40) 

k=l 

as the /c-point resolvent differential. As we can see, all are related to the fc-point resolvent functions except 
for W^"' 0. Define the Euler characteristics x of w^^^ aj^ 

X (u^r) = 2 - 2/1 - fc. (3.41) 

The sheet structure of this model makes it convenient to use the global loop insertion operator p.l7p to 
write the observables as differentials on the algebraic curve 



w 



'''\pu..:Pn) = (^ndx(p,)^j ^(''), (3.42) 



®Note that whenever the superindex in W^^"^ (x) distinguishes is as a term in topological expansion and thus is different from 
Wi{x) which is not an averaged quantity. 

^Noto that this is exactly the Euler characteristic of a surface with h holes and k boundaries. This allows us to make connection 
with the fatgraph interpretation of the topological expansion in which w^^'^ represents the generating function of counting numbers 
for fatgraphs of genus h with k boundaries. 
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where F^^'^ is the h order term in the topological expansion of F 



-^=E(¥) (3-43) 



h=l 

Finally define the differentials 



z.i'')^!, . . . ,p„) = ZJi'^bi, • ■ • + '5„,2<5.,OT^^v^^^fe^, (3.44) 

{x[pi) - x[p2)r 

which are just the same as w^^'^ above except for which is shifted by the double pole. 
3.5.1 Recursion equations 

If we collect the order h terms of the equations p.37p we obtain the following recursion equations for > 

i = l i = l 

3 



J2 viP^'^ )dx(p)«,f ' (p(^) ) = -iif") {x(p)) {Ax(p)f + i W^^^ {p^''> , p(^' ), (3.45) 

! = 1 i=l 



■ 1 '1 3 . , 

2 — 1 Z = l 1 = 1 



where 



h-l 



Wi''' {p, q) = 4'"'^ (P, 9) + E ^l"^ (P)^^'""^ 

rn— 1 

h — 2 h — m—1 

8.ndWi''\p,q,r)^wt'\p,q,r)+Y^ ^ ^y*™) (q)tZyl''— (r) (3.46) 

m— 1 n— 1 

/i-l 

. \ ^ / — (m) , s — (h — m-1)/ s — (™)/ N — (/i-m — 1)/ n — (™) / N — — m-l)/ n\ 

+ 2^ [^i {P)^2 {q,r)+w\ '{q)w)^ '[r,p) + w\ ' {r)w)2 {p,q)j- 

m— 1 

Notice that these are recursion equations indexed by the Euler characteristics x, in other words, the differentials 
u'^'''' are expressed in terms of other differentials with higher Euler characteristic x- 
We define as well the corresponding differential w|,''^ that will be used latter 

m— 1 

h-2 h-m-1 

Wt\p,q,r) = wi''-'\p,q,r)+Y^ w["'\p)w["\q)w[''-"^-\r) (3.47) 

ni=l n— 1 

h~l 

+ {'^i"'' iP)'>^2'~"^~^\'lyr) + wj"^(g)u;^''""~^^(r,p) + wj"^(r)u;^''~""^^ (p, q)^ . 
Finally we need to complement these equations with the fixed filling fraction condition at each order h 

<f^^^^4'\-) = 0- (3.48) 
As a corollary we obtain similar conditions for the n-point resolvent functions 

wit\;P2,...,Pn) = 0. (3.49) 

15 



A]' 



3.5.2 The two-point resolvent correlation function 

The "initial conditions" of the recursion relations above are the only data needed to solve the recursion relations. 
This data are the fundamental differential y{p)dx{p) of the spectral curve, the two-point resolvent differential 
''JJ2'\p: q) c^nd the filling fraction conditions for all h. The spectral curve is specified by the moduli, and we can 
obtain w^^p, q) by the action of the loop insertion operator on the fundamental differential y{p)dx{p), 

d ^ dilf\p)dx{q) 

^"^'^^ai^R^^^^^^"^^) = -^^ x{qy~^ ■ ^^-^^^ 

Notice that apart from a constant factor in front depending to which x-sheet p belongs, the sum in the above 
formula produce a second order pole in the diagonal x{p) = x(q). This fact, together with the vanishing A 
periods condition, indicates that 

dx{q)g^^-^^yip)dxip) = -^B{q^'\p) + ^ [B{q^"\p) + B{q'^^\p)) 

(1) , , 1 dx{q)dx{j>) ^^-^^^ 

where B{q,p) is the so called Bergman kernel. Remember that g^'^^ satisfy x{q^'''') = x{q) ~ x, and in this case 
we also precise g'*^' G Xk for convenience. In the last line of p.5ip we have used the relation 

tBiq^%p)= ';\'^'fl . (3.52) 
^ (x(q) - x{p)y 

We obtain similar results by applying dxiq)gy^^^^^ or ~dx{q) + 0V2(-x(q)) )- ^^^^^ operators 

combined in dx{q) g^^^) (see p.l7p ) give 



From this expression we find that 



4"'(g,p) ^ dx{q)^^w['\p) = dx{q)^^{U'{p) y{p)) 

^B{q,p) 



(3.54) 



{x{q) - x(j)))^ ■ 

Note that (|3.54p implies Ljj'^\q,p) = B{q,p). For latter purposes we write equation p.52p in the following form 



^'^\q,P^'y)^-J24'HqJ^^) 



(3.55) 



which can be trivially extended to 



wll'lj^{p^''\pi,--- ,p„) = -^wi'^i(p(^\pi,-- - ,p„). (-3 5g-) 

4 Solution of the recursion relations 

The equations (|3.45p obtained in the previous section are in fact non-homogeneous recursion equations. The 
non-homogeneous terms R^'^^ and D^'^^ are unknown meromorphic functions with at most a pole of order 1 
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and 2 respectively at x = 0. Although these functions are unknown they arc not an impediment to solve the 

equations. In order to take advantage of (j3.45p we will use the third kind Abclian differential dSp^oil) to rewrite 
{h) 

= V Resd5p,„(g)u;('')(p) ^^'^^ 

^ — ^ p— fa 

by using the bilinear Riemann identity, the condition p.48p and the properties of d5p.o(9)- The sum on the 
second line runs over all poles of the integrand except p = q, since the only poles apart from p = q are the 
branch-points on the algebraic curve, the set of points A is the set of branch-points of the algebraic curve. 
We then use (|3.45p to express locally (p) in terms of known -by recursion hypothesis- differentials plus the 
irrelevant terms and dC"). 

In our situation we will distinguish between two possible contributions, the ones coming from simple branch- 
points (branch-points joining only two sheets of the algebraic curve) and double branch-points joining all three 
sheets. 

4.1 Contribution from simple branch-points 

In this section we compute the residue 

RcsdSp.o{qWi'\p) (4.2) 

on a simple branch-point a. Around a simple branch-point the differentials on the algebraic curve can be split 
into two components w['^\p) = [wi^\p)]i + [wi^\p)]-i satisfying the following property 

[w^^\^\p))]x = y[w[''\p)]x. (4.3) 

Obviously on a simple branch-point i?^(p) ^ p so that we can drop the super- index. In the same way we can 
write y{p)dx{p) ~ [y{p)dx{p)]i + [y{p)dx{p)]^i. Suppose the point p is in the neighborhood of the branch-point 
a, it is clear that '0{p) is also in the neighborhood of a. The third point on top of x{p) (that we call p) is not 
in the neighborhood of a. The differentials w^'''' have poles only at branch-points. The terms depending only 
on p are holomorphic in x{p) when p ^ a. With all this we have 

y(j))dx{p) = [y{p)dx{p)]i + [y{p)dx{p)]-i , w[''^ (p) = [w[''\p)]i + [w[''^ , 

y{-&{p))dx{p) = [y{p)dx{p)]i - [y{p)dx{p)]^i , w^^''' {'d{p)) ^ [w[^\p)]i - [w[^\p)]^i , 

y{p)dx{p) = -2[y{p)]idx{p) , w^^'\p) = -2[w^^'\p)l , (4.4) 

3 „ 

^y(p('))dx(p)<)(p«) = ^y{p)dx{p)w[''\p) + 2[y{p)dx{p)U[w[''\p)U , 



i=l 

SO that in the neighborhood of the branch point a 
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Now we can introduce (14.51) into (|4.2 ''^°l 

Rcsd5,,„(g)<')(p) = Res ^^^^%^iL-^^y(pW)dx(p)z«('''(pW) 
P-i-Q p^a [y[p) — y[tf[p)))ax(p) ^ 

= Res W^'^^bW,p«) (4-6) 

= - Res ^^^pAi) W^''\p^ i){p)). 

p~^a [y{p) - y{-d{p)))dx{p) ^ 

The topological recursive expansion of the loop equation p.45p has been used in the second equality. Some 
terms, included have disappeared due to the residue operation and the assumption that there is no pole 

structure in y{p)dx{p) or w(^\p) when p — J> a. Notice also that in the last expression, W^2^ has been changed 
into Wi^^ by using (|336| . 

4.2 Contribution from a branch-point with rif, = 2 
In this section wc compute the residue 

'RssASp,o{q)w'('\p) (4.7) 

p— >a 

on branch-points with rib ^ 2. The contribution around these branch-points can be found similarly as with the 
simple branch-points. In this case however we have three components for every differential. 

y{p)Ax{p) = [y{p)Ax{p)]i + [y{p)dx{p)]0 + [y{p)Ax{p)]02 , w^l'\p) = [w^l'\p)]i + [wf'\p)]0 -f- [w''l'\p)]gi , 

3 3 

= Y.y{p^'^)dx{p) = Z[y{p)Ax{p)]^ , = = i[wt\p)]^, 

i=l i=l 

(4.8) 

where 9 = e~ is a root of unity. In order to find Wi^\p) = [w'^i'\p)]e + [w'^^\p)]g2 we present the equations 

3 



Y,v{P^'^)Mp)w^i\i'^) ^i{[y{p)dx{p)]o[wf\p)]e2 + [y{p)Ax{p)]e.[w^^\p)]e) 

i=l 

3 

5](y(p«)dx(p))2<)(pW) =my{p)dx{p)]ef[w^^\p)]e + ([y(p)da;(p)],2)2[^W(p)],.) 



(4.9) 



i=l 

and their solution 



[^^i\p)]b = ,1'^,, ,3 ( ([2/(p)dx(p)],.)2^2/(pW)dx(p)u;f)(p(')) 

([y(p)da;(p)]e2)3 - ([y(p)da;(p)]9)3 ^ ^ 

3 N 

-[y{p)Ax{p)]eY.{y{p^^)Ax{p)fwf\p'^'^) 



3 

-([2/(p)dx(p)](,)^^y(p«)dx(p)«;[''\pW) 



(4.10) 



-"^"Notice that [uij^'^' (p)] i as well as w^^\p) are holomorphic around the branch-point and do not contribute to the residue. 
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which gives the following representation for w[^'^ {p) around a ri^ = 2 branch-point 

1 f (b(p)dx(p)],.)^-([yb)dx-b)].)^^^^^(,;)^^^^^^^(.)^^(,)^ 



\i[y{p)dx{p)]o2r - i[y{p)dxip)]gy 



([2/(p)dx(p)]e2)3-([y(p)d.T(p)]9)3 ^^^^^ ' ^^^^ 1 (4.11) 



/b)dx(p)^y(p«)da;(p)<'b«) 

+ ^(y(p«)dx(p))2<)(pW)) . 



If we introduce equations p.45p and define uj{p,-d^ (p)) — {y{p) — y{-d^ {p)))dx{p), we get the following 

+ E (^J'(P) + y (P''^ )) MpM'"^ (p<"' , ) + E ^3"^ (P''' . P'"' . P'^' (4.12) 
{y(p)R<-''1 (p) + Pt'') (p)) (dx(p))3 / ' (p, (p)) '(p,i52(p)) ^(p,i?i(p).^^^(p)) '\ 

where again we used p.56p in the last equality. When we introduce it into (|4.7p . the residue contribution 
reduces to 

p-^a p^a Uj{p,§'-{p)) P^a U>(p,§^(j))) 

_ dSpA<i)M''\p^^Hp),^Hp)) 

p^a ui(j),'d^{p))uj{p,'d'^{p)) ' 
where again the term with i?''') and L'''"' disappear due to the lack of pole structure of the integrand. 

4.3 Residue formula for u;^'*^ 

We can split the set of branch-points as A = Ai + A2, where Ai is the set of simple branch-points and A2 
is the set of = 2 branch-points -in our case there is just one such branch-point-. Putting together all the 
contributions we have 

= -V V V Res ^SpAq)M"\p.^'{p)) _ dS,A<l)wt\p,^'{p)^^Hp)) ^^'^^^ 

^ ^ ^P^a Ojip.dHp)) ^ p-^a U!(p,-d^(p))uj(p,-d'^(p)) 



Note that the non-homogenous terms i?^'') and I?^'''' have disappeared under the residue operator. Also we 
recall the fact 
characteristic. 



recall the fact that those equations define all the w['^'' in terms of other such differentials with higher Euler 



4.4 Diagrammatic representation of the residue formula 

There is a very compact diagrammatic representation of the equation above. The diagrammatic representation 
is a very clear and simple way to write the equation (|4.14p , to understand the important features without getting 
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lost in the details. The reader not interested in the diagrammatic representation can skip this section. 
Define the diagrammatic vertices 



^ ^ ^ V^a Uji'p.'dH-p]) 
1=1 qGA, j=l 



< 



V Res 

' ^ n — ^rv 



dS'p,o('7) 



(4.15) 



where the black dots mark the leg carrying (p) and it is accompanied by an index when there is ambiguity. 
We define as well 




(4.16) 



After all these definitions the representation of equation (|4.14p is 




(4.17) 



4.5 Residue formula for w 



Define the following multi-differentials 



i=l a£Ai j = l 



d5t.o(q)>Vf2(t,l5^"(t),PK) 



w(i,z9^-(0) 



V Res 



u{j>,d^{jp))u{j>,d^{p)) 



aGA2 



(4.18) 
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where the notation is Pk = {pi, ■ ■ ■ ,Pk} and W2 ^| and W3 ^ are expressed as 

JCK m=0 

wi;;>{t,t',t",PK)=»«f+-,'\t,t',t",PK)+ E E HTi{*.pj)4-7;2~'^{*'.*"-pk/j) 

JCK m^O 

I ij^) f<l T-» \ (^ — m— 1)/,// , \ I (fi-) /,// T-» \ (/i — m— 1)/, ,/ n \ 

+ ,Pj)K-j + 2 (* .*:Pk/j) ,Pj)u'fc_j+2 (*:*.Pk/j) 

/i — 2 h. — m — 1 



(4.19) 



E E E ^^i(*-p.J^:il(^^p..K'!:;-:il{t^PK/.,u..)^ 



Jl,J2CK?7i=0 71=0 

JinJ2=0 

with K = {1, . . . , fc} is the set of index values and J, Ji and J2 are subsets of j, ji or j2 index values taken 
from K respectively. Remember that the differential is defined to be equal to zero so that, for example, the 
double sum in the definitions of (t, t', Pk) does not contain the term J = 0, m = nor J = K, m = h. It 
is clear from the definitions above that w['^\p) = w['^\p). In fact, we expect wl!'\pi, . . . ,p„) = w\!^\pi, . . . ,pn) 
for every k. To prove it is sufficient to check that dx{p) gyi^p-^ w^k^ (Pk) = w^k+^iVf^^p). 

4.5.1 The action of the loop insertion operator on the cubic vertex 

Consider the action of the loop insertion operator on the cubic vertex ^ with a generic /(i,i?(*)) attached 

to it. 

El oil, ,t!V^''° (*)) ^^(P) '^"^ {u{t,-»^{t)))^ '-"^^^ dv(p) 

+ Res4^d.(p)Mi^^ (4.20) 

-± E Efn-f^SlyE E ±^^%:^^,{Biu..p).s..^^^^ 

i = laeA^j = l\ ^' '•"fc = l/3GAj^Z = l ' ^ " 

- Re. d5,o(,)/(t^ ^-(.)) _ jiS^ ant^ 

We will drop from now on the last term involving the action of the loop insertion operator onto f{a,b). We 
take care of it by chain rule, after droping this term we can bring last equality to a more meaningful form if 
we exchange the order of the residues on the first term of the third line. The only difficulty to commute the 
residues appear whenever a = /? e Aj , in that case we need to apply the relation 

3 

Res Res = Res Res - Res Res . (4-21) 

k=0 ^ ' 

Two situations appear: when a € Ai the contribution from the last term of (|4.2ip cancels the last term of 
(|4.20p . when a € A2 the cancellation happens as well but some contribution remains giving rise to a term of 
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the type 



The result of exchanging the order of the residues i 



i0 



dx{p) 



dV{p) 



Uj{t, 1?J{t)) ' 



i — 1 aGA^ J 



■ a;(t,i?J'(t)) 



^—f t^a uj(^t,'t)'-{t))uj{t,'tt''(t)) 



Graphically, equation (|4.22p is represented as 



(4.22) 



dx{p) 



dVip) 




+ 



(4.23) 

where the differential w^\q,p) = B{q,p) is represented by the non-arrowed propagator P Q . 

4.5.2 The action of the loop insertion operator on the quartic vertex 

Consider now the action of the loop insertion operator on the quartic vertex 



j2 ( Res ./y:f^:fSL MP) + Reg, ,^ J'jc^^ mp) 

Qg A2 



dV{p) t^c, uj{t,^'^{t))Lj{t,S2{t)) 



Res (.it,,\t))Mpf''^yy +.(t,.^w)d.(p) 



av{p) 



av{p) 



dV{p) 



y Res 



f{t,i)\t),i)\t)) ^ E . {b{u,p)AS,A^\u)) + i?(^'(«),p)d5.,„(.) 



d5't,o{g)/(t,i?i{t),i92(t)) 



Res 

t^a (i.j(t,,3i{t))a;{t,i32(t)))2 



{u^it, 1 p) - Bi^-" {t),p)) + a;(t, i?^ (t)){B{t, p) ~ B(^^ {t),p))) 



Res __ii^^iiM__d.(p).^^(*''''«''''«) 



(4.24) 

As with the previous section we will drop the term containing the action of the loop insertion o perato r on 
f{a,b,c) and take care of it latter. After droping the term proceed by applying again the relation (j4.2ip and 

commute the residues in the first line of the last equality, then the residues of the type Rest_j.ct X]i=i ^^^u^€^{t) 



^^In appendix [D] there is a detailed account of this calculation. 
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cancel exactly the last line of (|4.24p and we obtain 



dx(p) 



- y Res 



dSu.oig) 



aeA2 

2 k 



y y y-Res 



B(u,p) - Y: Res 



+B{-&\u),p) i - y Res 



d5t,o(«) 



The graphic representation of this equation is 



dx{p) 



dV{p) 






(4.25) 



(4.26) 



4.5.3 The loop insertion operator on equation (|4.14p 

' dV{p) 



If we apply dx{p) qyi„\ to equation ()4.19p . and use equations ()4.22p and (|4.25p we have 



dxij}) 



a 



dV(p) 



dS„,o(p) 



i=l qGA; i = l 



W^\'(t,i>^(t),PK) i3(.5'(u),(p)) 



} 



y- dg..,o(p) 



V V )r Res ^ 

fc=i i3eA^ i=i *^ >■ 



+S()?i(u),p)W^|;' «, Pk) + p)>V^|;' (u, Pk)) 



+d:r (p) Vy3<^'^ (t , ^ ^ (p) , (p) , Pk ) } . 



(4.27) 



This equation simplifies considerably by regrouping some terms 

2 A: 



dx{p) ^^^^ ^ w^^^ 



(PK) = -g^£gRes_g^^{<^)^(„,PK)«,f(.'(«),p) 

+da;(p) -|— Wf ' («, 1?' («), Pk) 



P) 



a;(n, i3i(u))ix)(M, i9''(n)) I ^ -^A ^ 



;3eA2 



(tf2(„))^W(„^^l(„)^Pj^) + dx(p)^^W(';',(«,tfl(«),tf2(„)) 



(4.28) 



which is of the form ()4.18p . One can check that 



y^'i^l+M-a\n),V^,p) =w^^^^{u,Vy^)wf\i)\u),p) + w^^^^{^\u),V^)wf\u,p) + dx{p)^^^ 
^^^^l+^{u,-d\u),-d\u),Vvi,p) =«;<''\«)w(^^'(i?i{«),i?2(„),PK)+u,f'(^i(«))W^^^'(i?2(„)^„^Pj_,) 



(4.29) 
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which finally proves that w'j^]^^{F-K,p) = dx(j)) gy(^p) w^"'' (Pk)- Consequently we prove that w^"''(Pk 



w^I'^^IPk)- Thus, the residue formula for w^''^(Pk) is 



(/!.) 



V Res 



aeA2 



(4.30) 



4.5.4 Diagrammatic representation of u;^''^ 

The equation ()4.30p has the following graphical representation 





EE 

JCK m=0 



h-1 



EE 

JCK m=0 




/i.-2 h-m-1 

E E E " 

Ji,J2CKm=0 ri=0 




(4.31) 

Using the diagrammatic equation we can see how the /i-order /c-point resolvent differentials w^.'''' combine to 
construct other such differentials. We can interpret that the residue in the vertices connect two or three free 
legs (either of the same or of different ones). In the process some handles may be created. 

We can also see that they combine in a hierarchical fashion such that the equation (|4.30p form a recursion 
relation in the Euler characteristic x- w^'' and w^'' are the differentials with highest Euler characteristic, they 
are -as announced earlier- the initial conditions for the recursion. 

4.6 Variation of w^!^^ with respect to the moduh 

Now the we have found the recursion formula (|4.30p for w^''^ we can study its variation with respect to the 
moduli of the problem ta- We expect the variation of w^''' with respect to ta to be 



(4.32) 
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We prove this equation by recursion. From the Rauch variational formula we have 

2 i 

^it.tl,.tt*-"-^(*'''^W) (4.33) 
=-^-\J2 E Efi^'^ T'°Jfl^ {B{t,-)B{^Ht),p) + B{^Ht),-)B{t,p)) 

= -Ja(W^°Hq,P,-)), 

which proves (|4.32p for fc = 2 and /i = 0. This is the first equation in the recursion. To prove equation (|4.32p 
for all k and h we need some more results first. We already have 

dt^B{q,p) = -Ja (dx{-)-^B{q,p)] . (4.34) 



dV{-y 



Integrating this equation we have 



We also need 



dtASuAq) = -Ja ( dx{-)^^dSuA<l) ] ■ (4.35) 



B{t,-)-B{d^{t),-) 



--Ja 



Ja dx{-) 



d 1 



(4.36) 



'dV{-)u{t,d^(t)) 
We have now enough information to apply dt^ to equation (|4.30l) 

Pk) = g E ^(^^7(I» E^ 

(4.37) 

By recursion hypothesis we assume that 

dtM''\^^) = -Ja (M^i.+\(PK, •)) , (4.38) 
where x {^k^ > 1 — 2/i — fc. With that assumption we prove that 



dt^W^i;l{a, b, Pk) = - X ( dx(-)^wS(«' ^' Pk) ) , 

(4.39) 



dV{-) 

a*„>vg(a, 6, c, Pk) =-Ja (^d.T(.)^wS(a,^c,PK)) • 
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Using this result and (|4.34p . (|4.35p and (|4.36p into equation (|4.32p we have 



/_ ^ d5t.„(g)W^;'2(f,tf^(f),^;)^(t),PK) ' 
^ ^ ^^(•) V ctl.'^^" ^P,^i(p)Mp,'52(p)) // (4-40) 

= -Ja(wi%{q,PK,-)), 

which proves equation (|4.32p . 

Remark: Notice that in all these calculations the integral operators commute with the residues at the branch-points. This is 
normal since any of the contours in the definitions of Ja never approach the branch-points. 

5 Topological expansion of the free energy 

The topological expansion of the free energy = X^ft^o iw)'^'^ ^® related with the /c-point resolvent differ- 
entials w^!^'' by equation p.42p . Consider for the moment the simplest case fc = 1 

<'(p) = -dx(p)^^W. (5.1) 

By inverting the loop insertion operator we can find the coefficients J^^'*' of the free energy topological expansion. 
In this section we want to prove that there is an operator that we call Hx(p) such that 



w[^l,{V^,p)\={2-2h-k)wi^\V^) 



for all w''^^ with x < 0. We want to prove as well that 



w'('\p)\ for/i>l, (5.3) 



and so that the -ffx(p) operator is indeed related to the inverse of the loop insertion operator on all resolvent 
differentials with x < 0. 

5.1 The Hx(^p) operator 

Consider the operators defined in p.34p and construct from them the _ffj.(p) -operator 

H^{p)Wp)) = - X! ^'^^'^(V'(P)), (5.4) 

aeM 

where ipip) is a- meromorphic differential on the algebraic curve with pole structure only at the branch-points 
without residue. The Hx[p) operator is an essential piece of the whole construction of the topological expansion 
of the free energy, and as such has appeared very often in the literature [71 [TOl |8] . It has been proved in the 
literature that the H^(^pj operator can be written as well as 

H.(P0(P)) = - V Res ^ip)i^{p), (5.5) 
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where a are all the poles of and ^'(p) is a primitive of j/(p)dx(p) 



y(p)dx(p). 



(5.6) 



In equation p.35p we saw that Hx(^p^(B{p, q)) ~ —y{q)dx{q). For our purposes this result is better expressed as 



H, 



x{p) 



w^2^ip,Q) -y{q)<ix{q). 



(5.7) 



This is not of the form (|5.2p because w^"^ has a positive x- Notice that i?a;(p) is an integral operator that acts 
only on the p- variable. 

We use the equation (|4.30p to prove (|5.2p by recurrence. First we prove that it is true for the first term of 
the recursion, in our case this term will be Next we prove (|5.2p for generic (fc, h) assuming it to be true 

for resolvent differentials with higher x- 



5.2 The Hx(p) operator acting on wf^ 

The first equation of the recursion (|4.30p is for w^^^ wc have to apply Hx(^p^ to that equation to check equation 
(lOl) on w^°^ 



H. 



x{p) 



wf\pi.,P2-,q) 



--H 



x(p) 



i=l aeAi 3 = 1 ^ ' ^ 



1=1 Q6Ai j = l 
2 i 



(\.St,o{pi) 



{B{t,p2)y{dHt))dx{t) + B{d^t),p2)y{t)dx{t)) 



E E?-i^^B(.^w,p.).(Od.(„ 

1=1 a£Ai j = l ^ ' ^ " 



(5.8) 

Each of the residues on the third line vanish due to the lack of pole structure of the integrand. Notice that this 
result satisfies (|5.2p because x (^3°'') ~ 0- 

5.3 The -ffx(p) operator acting on w^,'*^ for k > 2 and x < 

Consider with x = 1^2ft, — /c: 



ip) 



w[%{q,Pli,p) 



y y y Res 



i=l aSAi j=l 



y Res 



dSt,,iq)W^,''^{t,^\t),i}\t),Pi,,p) 



t^a Uj{p,^\p))uj{p,'d^{p)) 



(5.9) 
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Next we introduce equation (|4.19p to make the differentials w!^'' explici10 and express (|5.9p as the sum of three 
terms 



y V y Res 



(^(°) (t, p)^^';), it), Pk) + (^^^- {t),p)w[''^^ it, Pk)) 



ERes 



dSt,oiq) 



(^(P,1?1(P))'^(P.'?^(P)) 



+ J/, 



x(p) 



+wi°\&\t),p)wi';l{^Ht),t,F,^) + ^i°\^Ht),pM'!l(t,&\t),FK)) 



■E E ER-^!?^dx{p)^w(;;'(M^{t),PK) 



(5.10) 



1=1 qS A; j = l 



a:{p) 



. Res- -l%4^1__dx(p)— ^W^f '(i,,3i(t),il2(f),Pj^) 



Let us treat each one of those three terms separately. The first one needs to be further expanded using equation 

l+k\ 



6301) for w['2kici,PK) 



i=l oeAi 3 = 1 



dgf,o(g) ,(0) 
ui{t, i9J(t)) 



(t, p)t«<';', (i^^^ (t) , Pk) + ^-J"' (*) , P)«'J + fc (*, Pk )) 



^ Res 



dSt,„(<;) 



(.c;(p,i*i(p))"(P,'?^(p)) 



+ tO<'"(lj'(t),p)W<''2()?'(t),t,PK) +«4°'('''(*)^P)W^!l'(t.'''(*).PK))] 



E E EE- 



i = l neAj j = l 



i.i;(t,i53(t)) 



V V V Ros """'"^^ 



>V^\;(u, !?'(«), Pk) 



(5.11) 



ERcs 



dSt,„(9) 



(<^t,p)W<;','('?'(*)^ '?"(*). Pk) 



(i^i (t), P) W<;;) (i*" {t),t, Pk) + ' (<3' {t),p}wi':l it, it), Pk)) 



■E E EE- 



i = l qEA; j=l 



Lj(t, l?J (t)) 



3 w(M,l5l(M))w(tl,)?2(„)) 



W^;fc'(u,i;^(«),i;"(u).PK) 



+ <'(i;'(t),p) - E 



/3eA 

dS^,o(t) 
^5 w(u,i3i(«))w(ti,i92(„)) -3,fc 



W<\'(«,'''("),'5'(").Pk) 



By applying equation (|5.7p . commuting the residues and after some lengthy algebra (sec appendix [Dl for a more 
detailed report of this computation) we get the following contribution 



+ 2 E ( tfr/tr .. m,u{u,^\u),^\u),V^). 

^ ^'^p uj{u,v'-[u),'iy^{u)) 



(5.12) 



The second and third terms of (|5.10p are much simpler. Notice that by the recursion hypothesis we assume 

Ah') 



53) to be true for x{wl^ ^ > (1 - 2/i - fc). It can be seen that all the resolvent differentials in the expression 

dVip) ' 



dx{p) gy, -. W^^fc* have an index x ^ r — 2h — k for ?' = 2, 3 and so we can use (|5.2p for them. Applying repeatedly 



-"^^Note that the only fe-point resolvent differential in W^'j,' that do not satisfy 115.21 1 is Making them explicit we can treat 

them separately. 
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(15.21) we have 



H. 



d 

x{p) ^dx{p)-^^^wi'^l{a,b,PK) 
da^WlJTTTrrWf^K^: c,Pk) 



= -{2h + k-2)W^%{a,b,PK), 
= -{2h + k-?,)wll{a,b,c,V,i), 



ip) 



dVip) 



Collecting all the contributions from (|5.12p and (|5.13p back into (|5.10p we have 



x{p) 



(h) 



{q,PK,p)] =-{2h+ik + l)-2)i -Y^T. El^es 



= -(2/i+(fc + l)-2)4';\(g,PK), 
which by recursion proves equation (|5.2I) . 

5.4 The operator acting on w[^^ for x < 0- Residue formula for T^'^^ 



Our goal in this section is to prove equation (|5.3p . 

Consider the following hierarchy of constants (depending only on the moduli of the matrix model) 



1 



2/1-2 



-H, 



x(p) 



-1 



y Res*(p)u;f'^(p) 



2/l - 2 ^ P^a 



(5.13) 



(5.14) 



(5.15) 



for ft, > 1, where as usual a are the branch-points of the algebraic curve. We have used here the form (|5.5p of 
the i?a;(p) operator. The claim is that these constants T'--^^ = J^^'^^ are equal to the coefficients of the topological 
expansion of the free energy. As a first check we act with the loop insertion operator on J^'''' 

{2h~2)dx(p)-^F^^'^ = -dx{p)-^Y Res-^{q)w^l'\q) 



dVipY 



dVip) 



dV{p) 



dV{p) 



The action of the loop insertion operator on the first term is 



dx{p) 



d 



dV{p) 



d 

dx{p)^^^y{r)dx{r) 



3 {x(j)) — x{r)y 



-dSq^p) 



dx{p) 



and so the result of the residue is 



VRcs 



( 9 



x{p) — x{q) x{p) — x{o) 



^{q)]w['\q) 



q^a \dV(j)) 



4=0 



//'\P) 



(5.16) 



(5.17) 



(5.18) 
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after moving contours and integrating. The second term of (|5.16[) is 

^Res*(,)d..(;,)^z.('"(,) 



x(q) 



(5.19) 



d 



And collecting both contributions 



dx{p)- 



-w 



{h) 



dVix) 

This is telling us that ^'''^ — J^C*) may only depend on T and the filling fractions e. 
In fact we can work out any variation with respect to the parameters ta G M, 



(5.20) 



(2.-2)|-^('^) 



VRcs 



Ja{B{r, ■))w^^'\p) + ^ Res *(p)J, U''\p, ■) 

-Ja (^^Resd5p,o(-)«ii''^(p)^ -Ja (i/,(p)^f'(p,-)) 

(2/i-2)J, (wj").), 



(5.21) 



dta 



where we have used the relation -f^wj^ ~ — Ja yuTklij proved above. These are expected equations for T^^^i 
so we get to the conclusion that — J^C') docs not depend on any parameter of the model and so 



2/1 - 2 ^"(P) 



(5.22) 



5.5 Homogeneity of T^'^^ 

Now we can check the homogeneity equation for J^^''^ 



(p) 



w[''^ (p) 



(5.23) 



which implies Y.a tadt^T'-''^ = {2h - 2)J^('*) 
5.6 Commutation relation of Hx(p) and dx(g)- ^ 



'dV{q) 

In the space of x < resolvent multi-differentials these two operators satisfy 

d 



(p) 



dx{q) 



dViq) 



(5.24) 
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where I is the identity operator. Indeed, define Wq^^ ~ — J^^''^ and consider the commutator of the i?a;(p) operator 
and the loop insertion operator acting on a generic resolvent multi-differential w^f!2i with fc > and x < 0, 



ip) 



dx{q) 



d 



dViq) 



(h) 



(p) 



which proves the statement. 



dx{q) 



d 



dV{q) 



(h) 



i(p,Pk) 



dx{q) 



{2h +{k + l)- 2)u4';\(g, Pk) + (2/1 + k~ 2)wlX^{q, Pk) 



[h] 



(h) 



i(<Z,Pk), 



(5.25) 



5.7 Computation of leading and next to leading orders of the free energy. 

The leading order and the next to leading order of the topological expansion of the free energy cannot be found 
through formulas like (j5.3p . they have to be computed independently instead. Fortunately in our case not 
much is changed from the case already studied in [T]. Indeed, The leading term (as with all i(;^°^(Pk)) is 
not changed by the convergence of the two branch-points, and as for J-'^^\ the recursion formula for w[^\p) is 
identical to the usual Eynard-Orantin formula, and thus J^'^' remains unchanged except for the presence of a 
hard edge. For the sake of completeness we include them here without its derivation, which can be found in the 
literature (see [1]). 

5.7.1 Formula for J'^") 
Following [T] we have 

=Y,taJa{y{p)dx{p)) =Y.^atbJa®MB{p,q)). (5.26) 

where the operators J'a are defined by equation p.34p . The reader will find more information in [T] and [50] 

5.7.2 Formula for 

The first order contribution to the topological expansion of the free energy is almost the same as the one found 
inH] 

(5.27) 



where y~'^{(3) is the coefhcient of z{p)^'^ in the jet expansion of y{p) around /? £ A2I1I y{p) = y~'^{(3)z{p)^'^ + 
{(3)z{p)^^ +y'^ {l3)+0 {z{p)) where z{p) is a local coordinate around (3. The only difference with the expression 
in [1] is that the product over the branch-points does not include the hard edge at /3 whose contribution is given 
by the (j/~^(/3))* term. The derivation of this term follows exactly the same lines as that of the other terms. 
The contribution of the double branch point (3 to tu^^"* is 



where B{f3, q) 



B(p.q) 



P^P 



' 9 j/-2(/3) ' 

and z(p) is a local coordinate around /3. Looking at the jet expansion of the 



equation (j3.53p wc obtain at leading order the relation 



3da;((7) 



d 



dV{q) 



y-\[3) = -B{P,q)) 



(5.29) 



This proves the contribution {y ^(/3))^ in J'^^). 



^^It has been already pointed out elsewhere that y{x) oc x as x — > whenever x = becomes a hard edge for both Mi and 
M2 (see for example 1181 [T)), thus the jet expansion must begin at order z{p)~^ 
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6 Conclusion 



In this article the topological expansion for the free energy and all fc-point resolvent correlation functions of the 
Cauchy matrix model with one non-simple (but non-critical) branch-point is successfully derived. This model 
lies right outside the reach of the standard topological expansion where only simple branch-points arc allowed. 
The new and larger structure that is uncovered here contains the Eynard-Orantin topological expansion formulas 
and extends them to a more general situation. 

Though the results of this article involve only branch-points with rib = 2, the possibility of generalization 
to the situation where the multiple branch-points are of any order seems straightforward. Those equations are 
not written yet but an immediate guess can be formulated from the results here and will be precisely derived 
in a future work. 



A Formal Matrix Models 

In this section we present a simplified definition of Formal Matrix Models. For more complete definitions see 
Consider the matricial S{Mi, M2) action. The action of the Cauchy matrix model is 

S{Mi,M2) = ^tr {Vi{Mi)) + ^tr (1/2(^2)) + tr (ln(Mi ® I + I ® A'h)) . (A.l) 

The formal matrix model is defined as a controlled perturbative expansion of the action around one of its minima. 
For this purpose we define M* and as the value of the matrices on one of these minima. Thanks to the 
U{N) invariance of the action (and the measure) we can supose M* to be diagonal and to avoid complications 
that would only obscure the definition will present only the one-cut situation (or genus zero case) in which M* 
is proportional to the identity. We refer to the above cited works for a more complete definition. Around this 
particular minima we have 

S{Mi,M2) = S{M*,m;) + -^{SM^f -f —^{5M2f + 5S{MI,m;, 5Mi,SM2), (^-2) 

where under the assumptions above we have 6 Mi = Mi — m*I and m* is the value of the unique eigenvalue 
of M* and ^ is the coefficient corresponding to the quadratic term in the action. For simplicity we drop the 
dependence on M* of dS. Consider the finite gaussian integral, it is proven that such an integral is just a finite 
polynomial in jj 



dA/idA/2 {SS{Mi - mil, M2 - mil))" e" 



Note the combinatorial interpretation of this integral in terms of Feynman graphs gives meaning to the quantity 
as a means of grouping together terms with different topology, thus the name topological expansion. We call 
formal matrix model the sum 

^SiMtMi) y til: / dMidM2 [5S{M, - mil, M, - m*l))" ^-f (4^'(a/i-™Ii)^+4^'(a/.-™P)^) 



n 

n=0 



, (A,4) 
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where the n-sum has a finite number of terms for every vakic of I but the whole Z-sum is divergent. Note that, 
would we be allowed to exchange the integration and n-summation, wc would obtain the regular convergent 
matrix integral 

-S(Mi,M2 



dMidAf2C-^^^^i'^^^^ (A.5) 

As a means of simplifying notation we use expression (jA.5p instead of (jA.4p while explicitely stating that the 
nature of our research involves only formal matrix models. 

B Loop Equations and definitions 
B.l General Form 

The loop equations are a glorified version of the following fact: suppose we have a multivariate integral on a 
domain fi, f{x)dx and a vector field x ~ h{x). Then the action of the vector field on the integrand induces 
the integral of a "total derivative" , namely 

^ ^ .f{x)A^x ^ l^{f^-^f + /div(/J)) d^x = 1^ div (//j) = - f{x)h{x)d''-^x, (B.l) 

where the last integral is the "flux" of the vector on the boundary of f2. If either / or ft, vanish on (or are 
tangential to) the boundary, then the variation is zero and thus yields identities amongst different integrals. In 
our case x are the two matrices and the domain of integration Q is the cone of positive matrices. Consider the 
vector field induced by the infinitesimal change of variable 

Mj ^ Mj + eSMj . (B.2) 

The partition function Zc is an integral on the space of positive definite matrices; the infinitesimal variation of 
the integrand under (|B.2[) is the divergence of a vector 

(B.3) 

Zc^Zc + tbZc + O . 

In order to have a vanishing variation the vector field must vanish on the set of matrices with at least one 
zero-eigenvalue. Under this condition wc will have 8Zc ~ 0. From the explicit form of Zc ()2.1|) . SZc = is 
expressed as 

JiSM,)) - ^ (tr(a.fiF/(A/0)) - (tr. (sM.tr. L^iIi^mJ )) = ^^'^^ 



J{SMi) is the factor coming from the Jacobian of the change of variables. This is the trace of the matrix 

d(Mi)k,i ' 

5(Mi),,j 
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B.2 Quadratic and Cubic loop equations 

The equations p.4[) and p.7p were already derived in [1] . We are going to provide only the necessary definitions 
to read the equations and the infinitesimal change of variables originating the loop equation. We refer the reader 
to the cited article for details. 

For the quadratic loop equation (|3.4|) we need the changes of variables 



T / 1 1 



and dM2 



N \x- Ml 
T / 1 1 



N \-x~ A/2 -X 



and the corresponding Jacobians 



'T / 1 
J((5A/i) = ( — tr 



N \x- M 



T / 1 
and J{SM2) = ( — tr 



,7V \-x-M2 

and the definitions of the following objects (apart from Wi, Yi and Ui) 



N \ X- M 



N \ x + M2 
Ji = ^tr(y/(A/i)) + Str^ ^ 



N ^ " \Mi®I + I(E) M2 

l2^-lti:{V^{M2))~^t/ 



To get the cubic loop equation p.7p we need to consider 

^ N^\x-Mi X J \N \x-AIiJ N \-x-M2) ' 

SM2 = ^ f -] ( — tr ( ^ - — tr ( ^ - V'(^x) + V'(x) 

N^\-x-M2 -xJ\N \-x-M2) N \x-MiJ ' ' 

and the corresponding Jacobians 

- (S- (^))^ - (^))^ - - (S- 

T'i f 1 \ 1 / 1 

rtr — ^ - ^-tr 



3 / T' / 1 \\2 



T3 / 1 \ T^* 1 



Af3 V(r;-M2)^y Af=* »? V(r;-Af2) 



(B.6) 



(B.7) 



N ' " \Mi®l + l®M2, 

If we put all this together into (jB.4[) . after some simplifications we have 

{{Y,[x)f) + {{Y2{x)f) + {y,{x)Y2{x)) = (r{x)) , (B.8) 

where 

Rix) = ^{V{{xf + Vi{~xf - Vi{x)Vi{-x)) P,{x) P2{x) 
= {U[{x)f + {U'2{x)f + - P,{x) - P2(a;). 



(B.IO) 



(B.ll) 
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We also need to define the objects 



^ , , r2 /V{{x)-V{{Mi) f 1 



iV2 '\x-Mi \Mi®l + l®M2 

r. , ^ fyi{M2) ( 1 

S2{x) = ---tr2 -^tn ^ 



Putting all together we find three cubic equations, one for each Yi 

i^u{xf) - {R{x)nix)) - ^ + i (^w.ix)) = (Dix)) , for k = 0, 1,2, (B.12) 

where 

(Oix)) = - (Y,ix)'Y2{x) +Yiix)Y2ixf) 

= ~U[,{x)U[ix)U^{x) - U[[x) (P2(x)) - U^{x) (Pi(x)) - i [U[{x) l^h) + C/^ (/i)) (B.13) 
+ + Q2{x)) + (m^i(x)2m^2(x)) + (m/i(x)W^2(2;)') + {Si{x)) + (52(x)) . 

B.3 Pole Structure of D{x) 

In order to find explicitly the pole structure at x = we need a new loop equation not derived in [T]. Consider 
the following change of variables 



t3 / 1 1 MA / 1 
<5A/i =— ^ tr2 ' 



T3 / 1 1 M2 Y / 1 



iM2 =—T ^ tri 



\-x - M2 -X x"^ ) \M\®\ + 1® M2 

with the corresponding Jacobians 



(B.14) 



JiSM^) =— I tr I i ) tri I tr2 



+ -'12 t'l fi -.tr2 A/2tri 



A/i ® I + I ® A/2 / V A/i I + I ® Af2 J J ^ \ \ A/l ® I + I ® A/2 / \ A/i ® I + tt ® A/2 



,7(5A/2) =— tr ti2 trj 



, -a: - A/2 / V — A/2 \ ® I + I ® A/2 / / \ V —x - I ® A/2 A/^ ® I + I ® A/2 

+ — tri (tr2 ( ) tr2 ( ) ) " -^t^l ( A/ltr2 





1 




) ^""^ ( Ml g 


) II + II . 


® M2)) 



(B.15) 



V A/i ® I + I ® A/2 / \ M-i ® I + I ® A/2 J J \ \ A/i ® I + I ® A/2 / V A/i ® I + I ® A/2 , 

Introducing this change of variables into (jB.4p we get (after some algebra and lots of cancellations) 

~ Ic ((]V2*'^ )ti2 ( ^ J + J ^ j^j^ ) ) ) ^ ( N^*"^^ (^2(J^-f ( Ml (8 I + I ® M2 



\ ' ' ' V Ml I + 1 Af2 / / / \ "V " V Ml 8 1 + 1 «) A'/2 



{B.16) 



This loop equation specifies the pole structure of ( D{x) ) which amounts to a degree two pole at most 
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C Definitions on the algebraic curve 



In this appendix we will give a brief review of several definitions of objects on the algebraic curve. The notions 
are relatively standard and we refer for example to [HJ [52] . 

C.l The algebraic curve 

Given an algebraic curve E{x, y) = (irreducible polynomial equation of degree di in x and d2 in y), its points 
{x,y) parametrize the points p of a Riemann surface S. Choosing a; S C U {0} {y £ CD {0}) as the base 
coordinate there are ^2 (di) solutions yi{x) {xi{y)) of the algebraic curve. Those functions define ^2 x-sheets 
(di y-sheets) which form a chart on the Riemann surface. Since, for every point p G S there is a pair (x, y) 
that solve E{x,y) = we can define two mcromorphic functions x{p),y{p) on S. For a generic p — p^"-* G S 
there are d2 ~ I other points p(*\ i ^ I, . . . ,d2 — I such that x{p^'^'>) — x{p). Each one of these points lies in 
a different x-sheet. In general, on these points we have yij)'^''^) ^ yip''-''') for all pairs The same definitions 
can be made based on y-sheets but we do not use it in this work. 

C.2 The branch points, branch cuts and locally conjugated points 

We will consider curves where dx{p) and dy{p) never vanish at the same point p (i.e. without cusps). The points 
where dx{p) = are called branch points and are denoted by a^. We only consider algebraic curves where the 
zeros of dx{p) arc simple except at a; = where we have a double zero. On simple branch points, the sheets 
are connected pairwise by cuts that go from one branch point to another, every branch point having only one 
attached cut. The branch point at a: = joins the three sheets together and belongs to two different cuts. We 
call Ai the ensemble of simple branch-points and A2 the ensemble of double branch points In a neighborhood Ua 
of any branch-point a G A^, given a point p £ Ua there are i points . . . G S such that x{p^^^) = x{p) 
for j ~ 1, ... ,1. The map : Ua ^ Ua is the map defined moving the point x{p) £ C around the point 
x{a) = a a total i times ending at the same value x and looking at the corresponding trajectory on Ua- We will 
order the points p^*-* defined above as p*-*-* = [p) . Note that the operation i?^ is defined only local and has no 
intrinsic meaning at the global level. Obviously, for a G A^ we have the following properties: 



C.3 Cycles and genus 

In every Riemann surface we can choose a basis of cycles for the homology group. These cycles can be chosen 
to form a canonical basis, i.e. two types of cycles Ai,Bi independent which satisfy the intersection conditions 
independent which satisfy the intersection conditions 



A Riemann surface of genus g has exactly g pairs of conjugated cycles Ai,Bi. Usually one defines the homotopy 
group by choosing related cycles A and B. These cycles are equivalent to the previous ones with the constraint 
that all pass through a common point pq. Cutting E along A and B defines E which is isomorphic to an open 
region of the complex plane C. This is called the canonical dissection of E 



^^{§^{p)) = d^+^{p)^p' 



(^j-\-k mod i) 



(C.l) 



A • A = = B ., ■ B. 



(C.2) 
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C.4 Abelian differentials. 

There are three type of meromorphic differentials on an algebraic curve. The differentials of the first kind 
(holomorphic) form a vector space of dimension g and we denote them by dui. We can always choose a basis 
dui{p) such that it satisfy 

f du,{p) = 5ij. (C.3) 
With this normalization we define the Riemann period matrix as the matrix of B periods 

du,{p). (C.4) 



This matrix is symmetric and its imaginary part is positive definite 5(t) > 0. These differentials dui(p) are 
called normalized Abelian differentials of the first kind. 

The Abelian differentials of the second kind are defined as meromorphic differentials with poles but 
with zero residues. A basis for these differentials may be called dil^,'^'' (p) which is a meromorphic differential 
with a pole at p ~ q without residue normalized as follows 

dn['\p) ^ (-fcz,(p)-^-i + o(i))dz,(p), 



dn'^>ip) = o, 

in some local parameter Zq(p) such that Zq{q) = 0. 

Finally the differentials of the third kind are meromorphic differentials that have only poles of first 
order. Again we will choose a basis for these differentials consisting of differentials dSq^r{p) with only two poles 
of order one, say at p = q,r, with residue 1,-1 respectively, and with vanishing A cycles, i.e. normalized 
following 

dSqAp)-(^-^+0{l)^ dZq{p), 

dSqAp) - + O {1)) dzrip), (C.6) 



■ip) 



j dV(p) = o, 



Mi 

where the local coordinates are defined as above. The differentials in this basis is called normalized Abelian 
differentials of the third kind. 

C.5 The fundamental bi-differential and relations with the Abelian differentials 

We also introduce the fundamental bi-differential (Bergman kernel) B(j), q) which is a meromorphic bi-differential 
with poles only at p — q oi order two with zero residue, and normalized in the following way 

N dz(p)dz{q) ,„ , , 
{z(p) - z{q)Y 



j B{p,q) = for i = 1, 
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This bi-differential is fundamentally eonneeted with the Abelian differentials presented above. 

Bip,q)^B{q,p), 
d f ip) ^ Res Bip,q)f{q), 

B{p,-) = 2niduj{p), 
dni^\p) = ResB{p,q)zr{q)-'' = -ResB(p,g)z,(<7)-^ 

iBj Jr 
rt—q rt—s 



(C.8) 



J t—v J t—r 



I t—p J t—r 

In all these equations, the integrations paths that are not closed cycles are constrained to not cross the A 
and B cycles, in other words belongs to the cut algebraic curve along the cycles. 

C.6 Riemann Bi-linear indetity 

Consider rj and cj two mcromorphic differentials. Consider also fJ, one of the anti-derivatives of lo based on a 
point Po 

n= f uj. (C.9) 

Jpo 

Under these definitions the Riemann bilinear identity takes the form 



V Res 17(p)?7(p) ^ 7r-y2\ f ^'^P^ ~ f ^^P^ f 



vip) ■ (CIO) 



This expression is the key for (|4.ip . 

C.7 Rauch Variational Formula 

Changes on the algebraic curve induce changes on the objects we have defined so far. The Rauch variational 
formula expresses how these changes affect the Bergman Kernel. Consider a generic variation Sdn, 

^ " dn{t)B{t,p)B{t,q) 



SdnB{p,q) 



dy(t)dx(t) 

(c.ii) 

1=1 qGA, j=l ' ^ 

The situations in which we are interested are the variations with respect to the moduli of our model, in particular 
the with respect to the potential. In that case we have the variation can be written as 

^^^")M^^(^''') = "£ ^ i^^'a '"lt'nt)t [Bit,r)B{^^{t),p) + B{^^{t),r)Bit,p)] 

2 i (C-12) 

= -E E EE^g jfil [Bit,r)Bi^^{t),p)+Bi^^it),r)Bit,p)] . 
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It is now easy to find the variation of other objects like the third type Abelian differential dSa,b{p) 

dx{r)^^dS^ AQ) = - E E E l^es d'^i [B{t, r)dSuA^^{t)) + B{^^ {t),r)dSuAt)] ■ (CIS) 

^ ' i=l aGAi 3 = 1 ^ ' ^ " 

D Details of the computations on sections 14.51 and 15.31 



Some of the calculations of sections 14.51 and 15.31 are cumbersome. This appendix contains the details of these 
calculations. 

D.l Calculation on section 14.5.11 

The action of the loop insertion operator on the cubic vertex (see equation (|4.20p ) consists of two main terms 

dx(p)— ^^f-V y y Res 4^^iiM^/(t.i)^(t)) 

where we have overseen the action of the loop insertion operator on the function f(a,b). It is important to 
notice as well that the function /(a, b) is symmetric under the permutation of its variables by hypothesis. We 
need to commute the residues of the second line. When the both residues are located on the same branch-point 
we use the relation 

i 

Res Res = Res Res — Res > Res , (D.2) 

where a £ A^. Residues around different branch-points commute without the need of extra terms. When the 
branch point a belongs to Ai we obtain 

= ^'^^ —, JTTT \B{u,p) Res dSt,a{^(u)) + B{^{u),p) Res dSt,o{u) p. 3 

-g?" ^"l , [Biu,p)dSt.o{-d{u)) + B{i){u),p)dSt,o{u)] . 

we evaluate the residues in the last line and we find 

-?<='^ ^u' tm! ^^f, ^ '1'^"'°^''^, [B(«,p)dSt,„(^(«)) + i3(i;(«),p)dSt.o(u)] = 

t->a Uj{t,-d{t)) ii->t,iS(t) 

/(t,.?(t)) r dS^(t).„(g) ^ , d5t.„(g) ^^^^ ' 

^^^^ds,.„(,)/(t.^W) 

When the branch-point under consideration belongs to A2 the commutation becomes more involved 

i^F-"^ ( ^[I'ZtW ] E (Biu.p)dStMi>'H) + B{^'{u),p)dStA^^)) = 

dSt.„(tf' («))/(«, (t)) , ^ ^ dSt.„(».)/(t, tf^' (t)) ] 
Blii, p) > Res h -d(w (u), u) > Res (V) 4) 



'i2?--( ^[l'Zm\ )i2 f ",°l'u (Biu,p)dS,A^'H) + B{^'iu),p)dS,Au)) 

^Ti*^" V"(*^'''(*))/ ^ ^ 
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We evaluate the residues of the last line and simplify all the terms and we get 



dS„,o(g) 



— — Res 



^ Res 



tj(t,i9'(t)) 



w(t,i9i(«))'^(t.''^(*)) ^ w(t,)?i(t))w(t,i92(t)) 

dSt,^ „(<?)/(*, i?^(t))B('5^(t),p) 



Res 



- Res 

3 = 1 



Cj(t,1?l(t))w(t,l?2(t)) 



W(t,l93(t))2 

;'B(t,p)/(i;^(t),i*'(t)) + S(i*\t),p)/(i;'(t),t) + B(ij'(t),p)/(t,i>i(t))) 



Putting all the contributions together, equation (jP.ip becomes 



dx(p) 



av{p) 



i = l QEAi j=l 



w(t,)?J(t)) 



fc=i /SeA^, 1=1 



- E 

QeA2 



/3 i3'(tt)) 



skp)E E Eft 

i=i oeAi j=i 



i=l qGA; i = l 



dSt,o(«)/(t,1?^(t)) 



a;(t,iJi(t));^(t,i?2(t)) 



t.(t,tfJ(t)) 

B(t,p)/(tfi(t),i;^(t)) + i3(i;i(t),p)/(i>^(t),t) + B(i>^(t),p)/(t,i*i(t)) 



where some terms from (|D.4p and (jP.Sp have cancelled the last term of equation (|D.ip . 



(D.5) 



(D.6) 



D.2 Calculation of section [4.5.21 

The action of the loop insertion operator on the quartic vertex (see equation (|4.24p ) is composed of two main 
terms 



da:(p) 



d 



av{p) 



ERos 



^ ( g?g J(t!f(S^t^Ht)) ^ ^ ES|4^^(^KP)dS*,o(^'(«)) + B(^'(«),p)dS.,„(n)) (D.7) 



Res 



dSt,„(g)/(t,i?H«),^'(*)) 

(<^(t,l?l(t))tj(t,l92(i)))2 



(<^(t,i*\t))(B(t,p) - S(i*'(t),p)) +a;(t,i>'(t))(B(t,p) - B(-d\t),p))) 



where we do not consider the action of the loop insertion operator on the function is a symmetric function 
/(a, 6, c). We must commute the double residue on the second line. The only situation that deserves some 
attention is when both branch-points coincide 



Res /(«-'^'(')-'^'(«)) y-Res '^•^"•"^'^^ (B{u,p)dSt o(i>'-(u)) + B{'9\u),p)dSt o(u) \ = 
t^c<^(t,tfl(t))w(t,lJ2(t)) ^„^^ V ^ \ \ \ ,,P! t,o\ JJ 

dS„.o(q) 



ERos 



Biu.p) Res '^f''"^"'^"", + B(^'(n),p) Res """^-"^^^ 



Res-^4&:^^E 



w(t,iJi (*))"(*■ 

dS„,„(g) 



' — 1 (t) 



w(t,1?l (*))"(*, l'^(t)) 

B(«,p)dS,,o(i9'(«)) + S(i3'(«),p)dSt,„(«) 



(D.8) 
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If we evaluate the residues in the last line and simplify the expression we get 

~Sf'' n^^^^ ^2/^^^ Z-^ ^ ^77^ ( -B(«, p)dSt,„(i? («)) + -B(t5 (u), p)dSt,„(«) j 

u->iS~'(t) 



w(*,l?l(*))w(t,l?2(t)) \^ w(t,lS»(t)) 



w(t,i5'(*)) 



Res 



and using the relation 



l^(t,1?2(t))2 w(t,l5l(t))2 



w(i3i(t),i?2(t)) Vw(t,'?^(t))^ w(t,)?i(*))^ 

_ w(i5i(t),iJ^(t))(w(t,i?^(t)) + w(t,i5^(t))) 

W(t,l91(t))2w(t,l92(t))2 



we get the final expression 



(D.9) 



(D.IO) 



Res ■ 



u ou,,, u ,2m-;T, ^3 — ^^^^ (i3(«,p)dSt,„(i;'(«)) + B(,5'(«),p)dSt,o(«)) 

Note that this cancels exactly the last line of (|D.7p . Putting all contributions together the result is 



dx{p) 



d 



OV{p) 



w(t,l?l(t))l^(t,lS2(t)) 



dS„.o(g) 



= y^ Res 



dSt o(i?'(ti)) . / . dSt 

B(«,p)Res — h S ()?(«), p) Res 



' w(t,-iJi(t))Lj(t,iJ2(t)) 



(D.ll) 



(D.12) 



D.3 Calculation of section [5.31 

The right hand side of equation (|5.1ip is composed of two terms, we are going to treat each of them independently. 



D.3.1 First term 

The first term is 



-EE E-air (-E E ES^^«--<-'. 



dS„,o(t) 



y(t)dx(t) Yl E I^<'| 



(D.13) 



dSt,„(9) 



»I=A2 



aj(p, i5i(p))a;(p, i?2(p)) 



(y(t)dx(t)/(i?'(t),>?'(t)) + y(iJ'(t))dx(t)/(i^'(t),t) + y(ij'(t))da;(t)/(t,'?'(t))) . 



The next step is to commute the pairs of residues. When the two residues are evaluated on the same branch- 
point we use again formula ()4.2ip put into another form 



Res Res = Res Res + Res 7 Res 

fa u— fa It— fa i— fa it— fa ^ — ^ t^^^ (t) 



(D.14) 
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(D.17) 



Two situations occur depending on the branching number of the branch-point. When a S Ai we have 

= - ^HrS '^'r/t'f^ (.(*)d5„,„(^(t)) + .(.(t))d5„,„(*)) (D.15) 
u-^a t-^Q uj[t,x}{t)j 

/(m, )?(«)) dSt,o(g)da;(t) /,/,\n,c. 
~ ^'''^ ^ JTTT ^'^f, > 7-^7TTT— (y(*)dS„,o(i?(t)) + y(i?(t))dS„,„(t)) . 

u^a t— ^it,T9(it) Uj{t,V{t}} 

The first line gives zero due to the Rest^Q. Evaluating the residues on the second line gives 

fiuMuYjf dS^(„),„(g)dx(u) , dS^,a{q)dx{u) \ dSVo(9) 

When the branch-point has branching number ni, = 2 the calculation becomes more involved, 

ERes 'If'-y fy(t)d.(t) f-E Res ^E2L^^f^fiu,i>'iu))] +yi^^t))d.it) f-f: R- f /(u, ^' (u)) 

To commute the residues of the first line we use equation (|D.14p again 

-E R- ^^^irSE (ds„,„(.^ (*)),(*) + .(^^w)ds„,„(*)) . 

As before, the first line is zero due to the Rest^^Q. The second term, after evaluating the residues is 

- E R- ^^r^, E 'VtT (d5„,o(^^- (*)).(*) + .(^^- w)d5„,o(t)) 

= -E R'='' •^/"'!i'/"u2 (-dS,,i(„) j9)y('5'H) + dS„,„(g)y(i;'(«))) 

+ R°" ^ .l/fw ,2. ^^ (dS^i(„)o('?)2/('''("))/(", +dS^2,„)„(g)y('?'(«))/(«, 

= + y" Res d.S...o(9) _ Res dS„.o(g)dx(M) + 

y(u)da:(n)/(i?^(n),i?"(u))dS„,„(g) / 1 ^ 1 

+ Rcs , „, , , — -r— — ; „, , + ■ 



(D.19) 
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Notice that the last term cancel exactly the last term in (jD.13l) . Putting together contributions from the Ai 
and A2 branch-points we have 



Res 



D.3.2 Second term 

The second term is 



■ E E^ E - ^ (<' (- E^ Kes .^(.), .^(«)) 

+ <'(^^(t),p) f- ^ Res— — -/Kfl\u).^^«))^ 

=y y Res ^^^^ Lt)d.(t) f- y rcs — '^f-'"^^'^*", — /k^^h.^^m) 

Again, to commute the residues we just need to consider coinciding branch-points 

_VRes^^i-M^ [v(t) Rc= dS..o(^^(*))/K^^M.^^W) dS„,o(t)/(«. ^H-), '?^(«)) 



Rcs 



j(u,w (Uj,v (U)) ^ y(W'(u))ax(u) ^ M 4- dS ( )) = 



(D.20) 



(D.21) 



= -R;s , Efi_g- (.(*)dS„,o(^^(t)) + .(^^(*))dS„,o(t)) (D.22) 

/(n,,5i(«),i?^(«)) ^ dSt,„(g)d:r(f) ^ ^^^^^ r■,:>f,^^A^ 
- R5S — JiTTT E 7-^-777^ (!/(t)dS„,o(i5-'(t)) -I- y(i5-'(t))dS„,o(t)j . 

The first line is equal to zero once more and the second becomes 



dS„.„(g)/(«,^^M,^^(n))dr.(u) ^^^^^ 1 1 ^ (D.23) 



^ d5^.o(g)/(«,i?' ^ (.;(tt) - y(i?^(n)))da;(«) 

^2 Res d5,..o(g)/(".'^'("),^'(")) 
"-So w(m, iJi(«))w(u,i32(„)) 
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